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ABSTRACT 

We  have  detected  a  satellite  via  optical  interferometry  for  the  first  time,  using  a  16  m  baseline  of  the  Navy 
Prototype  Optical  Interferometer  (NPOI)  to  observe  the  geostationary  communications  satellite  DirecTV-9S 
during  the  ’’glint”  seasons  of  February-March  2008  and  2009  when  the  sun-satellite-NPOI  geometry  was  favorable 
for  causing  specular  reflections  from  geostationary  satellites.  We  used  the  US  Naval  Observatory  Flagstaff 
Station  1  m  telescope  to  generate  accurate  positions  for  steering  the  NPOI.  Stars  are  the  easiest  targets  for 
optical/infrared  interferometers  because  of  their  high  surface  brightness.  Low  surface  brightness  targets  are 
more  difficult:  if  they  are  small  enough  not  to  be  resolved  out  by  typical  baselines,  they  are  likely  to  be  too  faint 
to  produce  detectable  fringes  in  an  atmospheric  coherence  time.  The  16  m  NPOI  baseline,  the  shortest  available 
at  the  time  of  our  observations,  resolves  out  structures  larger  than  ~  1.5  m  at  the  geostationary  distance,  while 
a  typical  size  for  the  solar  panel  arrays  is  2  m  x  30  m.  Our  detection  indicates  that  a  small  fraction  of  the 
satellite  glinted,  not  surprising  given  that  the  solar  panels  are  not  accurately  flat.  Our  fringe  data  are  consistent 
with  a  two-component  image  consisting  of  a  1  to  1.3  m  higher  surface  brightness  component  and  a  significantly 
larger  lower  surface  brightness  component.  The  brightness  of  the  glints  (2U4  and  <  1U5  for  the  two  detections 
in  March  2009)  and  the  size  scale  suggest  that  the  compact  component  has  an  albedo  of  0.06  to  0.13,  while  the 
larger-scale  component  is  much  darker,  if  circular  geometry  is  assumed. 
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1.  THE  NPOI 

The  Navy  Prototype  Optical  Interferometer  (NPOI),1  a  collaborative  project  of  the  Naval  Research  Laboratory 
and  the  US  Naval  Observatory  in  partnership  with  Lowell  Observatory,  is  a  multi-element  long-baseline  optical 
interferometer  located  at  the  Lowell  Observatory  dark-sky  site  on  Anderson  Mesa,  near  Flagstaff,  Arizona  (see 
Fig.  1).  The  NPOI  was  designed  for  both  wide-angle  astrometry  with  milliarcsecond  precision  and  milliarcsecond 
imaging  of  stars,  stellar  surfaces,  and  circumstellar  material.  First  fringes  were  obtained  in  1994,  and  the  first 
six-way  interferometric  observations  were  done  in  2001.  An  aerial  view  of  the  NPOI  is  shown  in  Fig.  1. 

The  NPOI  is  a  group-delay  tracking  interferometer,  which  increases  its  on-fringe  efficiency  but  requires 
tracking  the  fringes  rather  than  scanning  through  the  delay  range  where  the  fringes  are  expected  to  be.  For 
imaging  observations,  up  to  six  array  elements  (siderostats)  can  be  used,  although  it  is  common  to  use  three  to 
five  elements.  The  imaging  array  consists  of  two  moveable  elements,  with  a  third  expected  in  the  summer  of  2010 
and  three  more  under  development.  The  astrometric  array  elements  can  also  be  used  in  imaging  observations. 
The  astrometric  array  consists  of  four  fixed  array  elements  equipped  with  laser  metrology  to  monitor  the  motions 
of  the  siderostats  as  they  track. 
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Figure  1.  The  Navy  Prototype  Optical  Interferometer  (NPOI).  The  north  arm  is  the  one  extending  to  the  left  of  the  center 
of  the  array.  Feed  pipes  lead  from  the  array  center  to  the  optics  laboratory,  which  houses  the  continuously-adjustable 
“fast  delay  lines.”  Extending  to  the  north  of  the  optics  laboratory  is  the  array  of  “long  delay  lines”  that  will  provide 
additional  delay  when  the  longest  baselines  are  commissioned.  Close  to  the  center  of  the  array  are  the  four  enclosures  for 
the  astrometric  array  elements,  which  are  fixed  in  position.  The  east  and  center  astrometric  array  elements  were  used  for 
these  observations,  along  with  one  of  the  moveable  imaging  array  elements,  E6,  visible  just  above  and  to  the  right  of  the 
east  astrometric  enclosure.  Photo  courtesy  of  Michael  Collier. 


2.  GEOSATS  AND  GLINTS 

Geostationary  satellites  (geosats)  have  become  indispensible  for  communications,  surveillance,  and  environmental 
and  weather  observations.  Direct  imaging  of  these  satellites  could  be  helpful  in  diagnosing  problems  when  they 
arise,  but  at  an  altitude  of  nearly  35800  km,  1  m  subtends  5  milliarcseconds  (mas),  so  even  the  ~  40  m  span  of 
the  largest  solar  panels  is  considerably  less  than  an  arc  second  in  size  and  is  at  or  beyond  the  resolution  limit  of 
all  but  the  largest  adaptive-optics  equipped  telescopes. 

A  few  of  the  many  geosat  configurations  that  have  been  placed  in  orbit  are  shown  in  Fig.  2.  Geosats  are 
faint — roughly  14th  magnitude — almost  all  the  time.  However,  some  types  can  glint  brightly  for  a  few  minutes 
a  day  during  twice-yearly  “glint  seasons”  lasting  about  a  week.  These  types  consist  of  a  central  bus  and  two 
extended  solar-panel  arrays.  The  solar  arrays  are  oriented  north-south,  with  the  normal  to  the  panels  pointing 
to  the  the  right  ascension  of  the  mean  Sun.  The  alignment  with  the  Sun  is  not  exact:  first,  the  Sun’s  true  right 
ascension  differs  from  that  of  the  mean  Sun  by  up  to  ±4°  depending  on  the  time  of  year,  and  second,  early  in  a 
satellite’s  life,  when  the  panels  are  most  efficient,  they  may  be  pointed  several  degrees  away  from  the  mean  Sun 
in  order  to  avoid  generating  excess  power  or  solar  heating. 

Satellites  of  this  type  can  glint  brightly  when  the  Sun  reflects  specularly  from  the  solar  panels;  it’s  as  if  we 
see  the  Sun  reflected  in  a  mirror,  although  the  panels  are  not  optically  flat.  Therefore,  these  glints  can  occur 
only  when  the  declination  6  of  the  satellite  as  seen  from  a  given  location  is  close  to  the  declination  of  the  Sun. 
From  Flagstaff,  geostationary  satellites  are  at  <5  =  — 5?6,  so  the  glint  seasons  are  somewhat  before  the  spring 
equinox  and  somewhat  after  the  fall  equinox,  or  approximately  26  February-4  March  and  7-14  October. 

During  glint  season,  a  satellite  with  panels  pointed  toward  the  Sun  may  glint  when  the  azimuth  of  the  Sun 
(which  is  below  the  horizon)  is  180°  from  that  of  the  satellite,  when  the  Earth’s  shadow  covers  the  satellite. 
Only  at  the  beginning  of  the  (northern  hemisphere)  spring  season  or  the  end  of  the  fall  season — when  the  Sun  is 
slightly  south  of  5  =  —  5?6,  so  the  Earth’s  shadow  passes  north  of  the  satellite — can  a  glint  take  place.  However, 
if  the  panels  are  pointed  away  from  the  Sun,  the  glint  takes  place  before  or  after  the  Earth’s  shadow  covers  the 
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Figure  2.  A  few  of  the  many  types  of  satellite  in  geostationary  orbit.  The  satellite  detected  in  this  study,  DirecTV-9S,  is 
similar  to  the  HS-601  type  shown  here,  in  that  it  has  two  long,  narrow  solar-panel  arrays.  Typically,  such  satellites  have 
the  long  dimension  of  panels  oriented  north-south,  with  the  normal  to  the  panels  pointed  at  or  close  to  the  Sun.  This 
figure  is  taken  from  Ref.  2. 


satellite.  This  is  the  case  for  DirecTV-9S,  the  satellite  that  we  detected  with  the  NPOI  in  March  2008  and  March 
2009. 

The  glint  light  curves  shown  in  Fig.  3  demonstrate  these  two  behaviors.  AMSC1  entered  the  Earth’s  shadow 
before  the  glint  reached  maximum  brightness,  while  DirecTV-lR  had  already  passed  through  the  shadow  when 
it  peaked. 


Figure  3.  Light  curves  of  two  typical  glinting  satellites.  The  gap  around  6  =  0°  is  due  to  the  satellites  passing  through  the 
Earth’s  shadow.  The  normal  to  the  solar  panels  does  not  point  directly  toward  the  Sun  for  either  satellite;  for  DirecTV-lR, 
the  tilt  away  from  the  Sun  is  large  enough  to  move  the  peak  of  the  glint  out  of  Earth  shadow. 

The  theoretical  peak  magnitude  mg(peak)  of  a  glint  can  be  calculated  from  the  apparent  magnitude  of  the 
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Sun,  the  area  of  the  solar  panels,  and  their  albedo.3  For  a  panel  area  of  100  m2  and  an  albedo  A  of  0.1,  mg(peak) 
is  about  —  3™4,  much  brighter  than  the  observed  peak  magnitudes.  The  primary  reason  for  this  relative  faintness 
is  probably  the  lack  of  local  or  large-scale  flatness  of  the  solar  panels. 

3.  INTERFEROMETRIC  OBSERVATIONS 

Observing  geosats  with  an  optical  interferometer  presents  us  with  several  challenges.  The  first  is  the  glints 
themselves.  Not  all  satellites  glint,  and  predicting  whether  a  given  satellite  will  glint  is  not  straightforward.  If 
a  satellite  does  glint,  the  glint  may  be  extinguished  by  the  Earth’s  shadow.  And  if  the  glint  happens  outside 
Earth’s  shadow,  the  time  during  which  it  is  bright  enough  to  observe  is  usually  only  a  few  minutes.  A  second 
challenge  is  the  size  of  the  satellite.  Size  scales  of  ~  1.5  m  or  more  will  be  resolved  out  on  our  15.9  m  baseline, 
producing  very  weak  to  no  fringes. 

The  short  duration  presents  a  third  challenge:  finding  the  fringes  before  the  glint  fades.  In  order  to  keep  the 
fringe  search  time  shorter  than  the  glint  duration,  we  must  set  the  delay  lines  within  ~  300  gm  of  the  correct 
position,  which  in  turn  requires  knowledge  of  the  satellite  position  to  better  than  one  arcsecond  for  the  15.9  m 
baseline.  The  published  orbital  elements  (the  “two-line  elements”  [TLEs],  available  at  http://celestrak.com)  are 
not  precise  enough,  and  the  satellite  operators  can,  and  sometimes  do,  move  the  satellite  to  a  different  longitude. 
In  addition,  the  orbits  are  not  perfectly  geostationary;  we  found  that  drift  rates  of  0'.'5  min-1  in  azimuth  and 
elevation  were  not  unusual. 

To  reduce  the  position  uncertainties,  and  to  monitor  glints  as  they  brightened,  we  did  our  own  near-real-time 
astrometry,4  using  the  US  Naval  Observatory  Flagstaff  Station  (NOFS)  1  m  Ritchey-Chretien  telescope,  about 
19  km  distant  from  the  NPOI  site.  In  each  of  the  fields  that  we  monitored,  we  used  the  TLEs  and  the  program 
The  Sky  (http://www.bisque.com)  to  identify  the  satellites  when  there  was  more  than  one  visible  in  the  20'  field 
of  view.  We  continued  observations  until  the  glint  became  too  bright  for  the  detector,  typically  20  min  before 
maximum. 

At  the  NPOI,  we  used  the  three  array  elements  at  stations  E6,  AE,  and  AC.  Because  of  the  beam  combiner 
configuration,  we  could  use  only  the  E6-AE  (15.9  m  length  at  93°  azimuth)  and  AE-AC  (18.9  m  at  —68°) 
baselines,  so  no  closure  phases  were  available.  The  complexity  of  coordinating  the  observations  at  the  NOFS  1  m 
telescope  and  the  NPOI  led  us  to  concentrate  on  a  small  number  of  targets,  typically  three  to  five  each  night. 

Our  first  observing  season,  in  2007  October,  produced  experience  and  adjustments  to  the  techniques,  but  no 
fringes.  During  our  second  season,  in  February-March  2008,  we  detected  weak  fringes  on  one  target,  DirecTV-9S, 
described  in  Ref.  5.  The  2008  October  session  produced  no  detections,  even  though  DirecTV-9S  was  one  of  the 
targets.  However,  in  our  fourth  season,  February-March  2009,  we  once  more  detected  fringes  in  observations 
of  DirecTV-9S.  The  signal-to-noise  ratio  in  these  2009  fringes  was  considerably  better  than  2008.  The  fringe 
power  spectrum,  shown  in  Fig.  4,  confirms  the  detection.  In  both  2008  and  2009,  we  detected  fringes  on  only 
the  15.9  m  baseline.  We  suspect  that  our  non-detection  on  the  18.9  m  baseline  stems  from  the  ~  25%  reduction 
in  the  squared  fringe  visibility,  V2,  that  results  from  the  its  greater  length. 

Our  detection  data  consist  of  fringe  amplitudes  in  15  wavelength  channels  on  only  one  baseline,  too  sparse  to 
make  an  image.  However,  the  trend  of  V2  with  wavelength  can  give  us  the  approximate  size  scale  of  the  glinting 
area,  although  having  data  on  only  an  east-west  baseline  means  that  we  have  size  information  only  along  that 
direction.  Fortunately  for  us,  geosats  of  this  general  type  are  oriented  with  the  long  dimension  of  the  solar  panels 
north-south;  if  it  were  east-west,  the  satellite  would  be  resolved  out  on  this  baseline. 

The  2008  data  indicated  a  size  scale  for  the  glinting  area  of  roughly  1.3  m,  but  only  the  six  reddest  (lowest 
spatial  resolution)  channels  had  V2  significantly  different  from  0.  The  2009  visibilities  are  significantly  non-zero 
in  the  reddest  11  channels  and  have  much  better  signal-to-noise  ratio.  They  are  nearly  constant  across  that 
span,  suggesting  an  unresolved  component,  but  the  amplitudes  are  small  (V2  ~  0.07),  implying  that  a  resolved 
component  with  a  significant  fraction  of  the  total  flux  is  diluting  the  fringes  from  the  unresolved  component. 

We  present  the  2009  data  in  Fig.  5  as  |Vj  rather  than  the  more  conventional  V2  in  order  to  show  the 
contributions  of  the  two  components  to  the  visibility.  (The  2008  data  are  also  consistent  with  these  models, 
but  the  lower  SNR  of  those  data  constrains  the  models  poorly.)  Since  we  do  not  have  phase  information,  the 
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Figure  4.  Fringe  power  spectra  on  the  15.9  m  baseline.  The  abscissa,  k,  is  the  number  of  fringes  per  delay  sweep  stroke 
length  measured  in  wavelengths.  For  the  2008  observations,  the  NPOI  stroke  on  this  baseline  was  two  waves,  while  it 
was  three  waves  in  the  2009  observations.  The  power  peaks  at  k  =  2  and  k  =  3  in  the  respective  data  demonstrate  that 
fringes  were  detected. 


relative  positions  of  these  components  is  unknown;  for  example,  they  could  be  a  patch  of  the  solar  panels  plus 
an  area  of  the  bus.  For  purposes  of  the  calculation,  we  assumed  that  the  areas  are  circular,  although  we  have  no 
information  in  the  north-south  direction. 
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Figure  5.  Observed  visibility  amplitudes  in  2009,  and  two-component  model  visibilities.  The  visibilities  \V\  (rather  than 
the  usual  V2)  for  the  two  detections  in  March  2009  are  very  similar,  and  can  be  fit  with  similar  models  consisting  of  one 
component  of  ~  1  to  1.3  m  (upper  dashed  curve  in  each  panel)  and  a  second  component  larger  than  ~3m  (lower  dashed 
curves).  The  full  curve  shows  the  total  visibility  of  the  model. 


Table  1  shows  the  parameters  of  the  model  fits  to  these  data,  as  well  as  the  peak  glint  magnitudes  in  the 
Cousins  R  band  as  observed  at  the  NOFS  1  m  telescope.  (The  detectors  at  the  1  m  and  at  NPOI  saturated 
at  the  peak  brightness  on  3  March.)  The  data  do  not  constrain  the  size  (L2)  of  the  larger  component  well,  so 
we  fit  the  size  of  the  smaller  component,  L\ ,  and  the  fractional  flux  F\  for  three  values  of  £2-  The  fits  from 
the  two  nights  are  quite  similar,  although  the  best  fits  for  3  March  have  a  slightly  smaller  compact  component 
that  contributes  a  slightly  smaller  fraction  of  the  total  flux.  For  both  nights,  we  obtained  slightly  better  fits  for 
L2  =  3.7  m  and  7  m  than  for  the  intermediate  value  of  5  m. 


4.  DISCUSSION 

Are  these  size  scales  consistent  with  the  glint  magnitudes  we  observed?  The  maximum  expected  brightness  of 
a  glint  can  be  estimated  from  the  area  and  the  albedo,  as  in  the  approach  taken  in  Ref.  3  (see  above).  That 
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Table  1.  Two-component  model  fits  to  March  2009  data 


Li  (m) 

L'2  (m) 

F1/F2 

xi 

Tor  (peak)  (mag) 

2  March  2009: 

1.3  7 

0.6  /  0.4 

1.3 

2.4 

1.2 

5 

0.45  /  0.55 

1.8 

1.4 

3.7 

0.7  /  0.3 

1.4 

3  March  2009: 

1.1  7 

0.45  /  0.55 

0.8 

<  1.5? 

0.6 

5 

0.25  /  0.75 

1.0 

1.2 

3.7 

0.55  /  0.45 

0.8 

approach  can  be 
magnitude,  to: 


inverted  to  estimate  the  albedo  A,  given  the  size  scale  L  of  the  glinting  area  and  the  glint 


A  = 


(1) 


where  too  is  the  Sun’s  apparent  magnitude  (— 27™ 2  in  the  R  filter),  6q  is  the  angular  diameter  of  the  Sun 
(0.009  rad),  H  is  the  distance  to  the  satellite  (about  35800  km  for  a  geosat),  and  L  is  the  diameter  of  the  glinting 
region,  which  we  will  assume  is  circular.  For  simplicity,  we  also  assume  that  the  two  components  are  spatially 
distinct  (i.e. ,  the  larger  component  is  not  a  halo  around  the  smaller  one). 

Vrba4  measured  a  peak  magnitude  to  =  2™4  on  the  night  of  2  March  2009.  The  next  night,  the  image 
saturated,  and  may  have  reached  to  =  lm.  The  two-parameter  model  described  above  fits  both  nights’  interfero¬ 
metric  data  about  equally  well,  with  similar  size  scales  and  a  similar  division  of  flux  between  the  two  components, 
despite  the  difference  in  peak  magnitude. 


For  the  2  March  glint,  we  infer  A  ~  0.06  and  <  0.004  for  the  compact  and  extended  components,  respectively. 
If  we  assume  that  the  peak  brightness  on  3  March  was  1™5,  the  corresponding  albedo  estimates  are  ~  0.13  and 
<  0.012.  For  comparison,  Vrba4  estimated  A  =  0.9  in  the  2  March  data  taken  at  the  NOFS  1  m  telescope 
and  pointed  out  that  for  new  solar  panels,  A  «  0.12.  It  should  be  kept  in  mind  that  the  albedos  inferred  from 
our  interferometric  data  assume  a  circular  geometry,  but  we  have  no  information  on  the  north-south  extent  of 
the  glinting  areas.  If  the  “compact”  component  is  extended  north-south  along  the  long  dimension  of  the  solar 
panels,  the  actual  albedo  is  smaller  than  we  estimate  here.  By  contrast,  the  albedo  of  the  extended  component  is 
quite  small.  We  speculate  that  shiny  multi-layer  insulation,  or  MLI,  covering  the  spacecraft  bus  and  scattering 
sunlight  in  many  directions,  may  be  the  cause  of  these  low  values. 


5.  NEXT  STEPS 

The  results  of  these  observations  suggest  several  directions  for  further  work.  Perhaps  the  first  is  toward  observing 
at  lower  spatial  frequencies,  since  our  results  show  that,  at  least  during  glints,  the  flux  is  mostly  from  extended 
structures.  The  NPOI  is  unlikely  to  commission  short  (<  10  m)  baselines  or  an  infrared  detector  in  the  near 
future,  so  using  aperture  masking  techniques  on  8  m  class  telescopes  may  be  the  most  direct  route  to  short- 
baseline  observations.  We  are  currently  exploring  this  option. 

Further  steps  include  interferometric  observations  at  infrared  wavelengths  with  larger  apertures,  as  well  as 
increasing  the  number  of  apertures.  In  the  infrared,  albedos  are  higher,  and  the  spatial  frequency  sampled  by  a 
given  baseline  is  lower.  Observing  at  these  wavelengths  with  larger  apertures  could  extend  the  observing  season 
beyond  the  current  glint  season. 

The  greatest  challenge  will  be  acquiring  data  over  a  wide  enough  range  of  spatial  frequencies  to  allow  imag¬ 
ing  of  these  complex  objects.  Detecting  weak  fringes  at  the  higher  spatial  frequencies  that  characterize  detailed 
structure  will  require  multi-element  interferometers  with  some  means  to  phase  the  array,  such  as  baseline  boot¬ 
strapping. 
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